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Abstract 

A study is made of an anisotropic Potts model in three dimensions where the coupHng depends 
on both the Potts state on each site but also the direction of the bond between them using both 
analytical and numerical methods. The phase diagram is mapped out for all values of the exchange 
interactions. Six distinct phases are identified. Monte Carlo simulations have been used to obtain 
the order parameter and the values for the energy and entropy in the ground state and also the 
transition temperatures. Excellent agreement is found between the simulated and analytic results. 
We find one region where there are two phase transitions with the lines meeting in a triple point. 
The orbital ordering that occurs in LaMnO^ occurs as one of the ordered phases. 

PACS numbers: 05.10.Ln,75.30.Kz,05.70.Fh,64.60.Cn 
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I. INTRODUCTION 

There have been many studies of the phase transitions of the Potts modeli (for a general 
review see Wu 1982^) , largely due to the richness of its physical content and its relevance in 
real physical systems^. While a large body of exact and rigorous results are now known, a 
number of problems particularly those associated with models with antiferromagnetic and 
multi-site interaction are still being investigated. 

The ferromagnetic, FM, case has been well studiecii*^^ and it is now believed that 
the g-state FM Potts model in 3d for q > 3 exhibits a first-order transitioii^. In systems 
which order ferromagnetically, it is known that the critical behaviour of the system near the 
critical temperature Tq is not affected by the nature of the underlaying lattice. The critical 
behaviour depends only on the dimensionality d and number of components of the order 
parameters. 

Antiferromagnetic, AF, Potts models have been shown to possess interesting and unusual 
properties. The ground state entropy is nonzero whenever the number of spin states is g > 2. 
The g = 3 model on square lattice has critical points only at zero temperature^'^'^iiSiii. In 
three dimensions, the evidence indicates the existence of a phase transition for g = 3, 4 and 
5 that is believed to be weakly first order, although the nature of these transitions has been 
uncertai n^^i^^i^^i^^i^^i^''' and a finite entropy^^ is found at T = 0. The effects of frustration 
in the Potts model have been studied by considering competing nearest and next nearest 
neighbour interactions ™. 

A three state Potts model in three dimensions allows for the possibility that the Potts 
states and the space coordinates are coupled. This arises physically where there is strong 
Jahn- Teller coupling of an electronic doublet typically from d-electrons coupled to the two 
dimensional lattice distortions with strong unharmonic terms as discussed by Kanamori^S. 
This leads to three orbits which are used as the three states of the Potts model. 



\x) = (3x^ — r'^)f{r), 




\y) = {^y'-r')fir), 




\z) = {3z'-r^)f{r), 


(1) 



where /(r) is a radial function and = + + z^. 

In such a model the interaction between orbits depends on both the type of orbit and 
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FIG. 1: The orbital ordering in LaMnO^. It is antiferromagnetic in x-y plane and ferromagnetic 
along the z direction. 

the direction of the bond between them. The compound LaMnO^ has orbital order of this 
type, as seen in Fig. Q and the interactions between the orbits have been calculate d^^i^^i^^ . 

We define the anisotropic three dimensional Potts model in terms of two interactions, Ji 
and J2 indicated in Fig. |21 The interaction Ji is for two sites occupied by the same orbit, 
say X where the bond between them is along x direction (a head to head configuration), y 
orbits in y direction or z orbits in z direction as shown in Fig. |2K- The interaction J2 is 
when the orbits, say x, are separated by a bond in one of the other directions (a side to side 
configuration) y direction or z direction, y separated by a bond in x direction or z direction 
or z separated by a bond in x direction or y direction as shown Fig. ^p. Finally, any two 
orthogonal orbits have zero exchange interaction as shown Fig. |2t- 

In this paper we investigate the phases over the whole J1-J2 plane. We find an extraor- 
dinarily rich phase diagram. There are no less than six distinct phases and special critical 
properties on the lines separating them. In the case where Ji = J2 = J we recover the 
results of the isotropic three dimensional ferromagnetic Potts models (for J > 0) and the 
antiferromagnetic Potts models (for J < 0). Also, the ordering of LaMnO^ occurs as one 
of the phases. We find that all of the phases except the ferromagnetic and isotropic antifer- 
romagnetic phase are frustrated as the Potts states on all the sites cannot be arranged to 
optimize all the interactions. 

We use Monte Carlo simulations to identify the ordering that occurs in the low tem- 
perature limit in each of the phases and find the energy and entropy of the ground state 
of each phase. The simulated values for the specific heat are used to find the variation of 



3 



(b) J2 




FIG. 2: The types of the orbital interaction used in our simulation (a) Ji refers to head-to-head 
ordering in one direction, (b) J2 effect is to order the same states in parallel to form layers in two 
directions, (c) orthogonal ordering is considered as a zero interaction. 

the transition temperature over the J1-J2 plane and to investigate quaUtatively the order of 
the transition by comparing the form of the specific heat anomahes observed with the well- 
studied cases of the ferromagnetic and antiferromagnetic Potts model. The methodology is 
described in section 2. In section 3 the results are presented for the six phases. Results for 
the phase boundaries are given in section 4. Finally the conclusions are given in section 5. 

II. METHODOLOGY 
A. The Model 

The Hamiltonian for the standard (isotropic) Potts model Eq. ^ and the anisotropic 
three state Potts models Eq. Q on a simple cubic lattice are given below. The factor of 
1/2 is included to correct for double counting. 

7 ^ 

(id) 

where Si = x,y ot z is one of the three states on site i, Ss^,Sj is the Kronecker delta which 
equals 1 when the states on sites i and j are identical, Si = Sj, and is zero otherwise, 
means that the sum is over the nearest neighbour pairs, J is the integral exchange and N 
is the total number of the sites in the lattice. The anisotropic Potts model, that is studied 
here, differs from the standard Potts model with Hamiltonian given by Eq. Q because the 



4 



exchange interaction depends both on the orbit and on the direction of the bond p, 

1 ^ 

Hais = " 2 5Z ^Pij )^s,,s,, (3) 

where p. . = Ri — Rj- 

There are two exchange interactions for the anisotropic model. The 'head to head' inter- 
action Ji is defined by Jx{,p) = Ji for p = ±xa, Jy{p) = Ji for p = ±ya and Jz{p) = Ji for 
p = ±za. The 'side to side' interaction J2 is defined by Jxip) = -h for p = ±ya or p = ±za, 
Jy{p) = J2 for p = ±xa or p = ±za and Jz{p) = J2 for p = ±xa or p = ±ya. Thus each site 
has a couphng J2 to four neighbours and a couphng Ji to two neighbours. This is shown in 
Fig. 121 It is worth mentioning that these types of interaction do not affect the overall cubic 
symmetry of the lattice. 

For each phase we follow the procedures below. 

1. We perform a Monte Carlo simulation to find the nature of the ground state order. A 
phase is defined as a region that exhibits the same configuration in the ground state. 
We identify the nature of the long range order and the fraction of the sites that form 
the ordered state. 

2. From the observed ground state configuration we obtain an analytic expression for the 
ground state energy per site Un{0) and for the ground state entropy per site s„(0) if 
possible, where n is the phase number. The energy, Am„, which is the ordering energy 
of the phase is obtained as Aun = Mn(oo) — w(0) for each phase. 

3. We compare the analytic values of the ground state energy and entropy deduced from 
the observed order with those from the Monte Carlo simulations for each phase. 

4. We confirm that the phase line between two regional phases with different ground 
state order occurs for values of Ji and J2 such that the ground state energies of the 
two phases are equal. 

5. The Monte Carlo simulations are used to estimate the transition temperatures, Tc, 
around the phase digram and the nature of phase transition (first or second order) is 
determined in some cases. 
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6. We use a combination of the Monte Carlo results and analytic results to obtain in- 
formation about the type of order and ground state entropy, sl„(0), occurring on the 
boundary lines. 



B. Monte Carlo Simulations 

Our Monte Carlo calculations have been carried out on 3d finite cubic lattices (with 
linear size L = 8) with periodic boundary conditions. All our simulations have made use 
of the Metropolis algorithm with the spin being chosen at random, and with averaging 
performed over 10^ Monte Carlo steps per site. In most of the phase diagram this gave clear 
results. Where convergence was slow we checked that we had found the true ground state 
by both increasing the number of Monte Carlo steps and also by looking at a L = 10 lattice 
as explained below. Results at low temperatures were obtained by cooling down from a 
high-temperature random configuration as discussed by Banavar et al^. 

The internal energy per site of the system obtained from the simulation is given as follows, 

uiT) = ^{Hap)t, (4) 

where = L^. We checked that u{oo) = — + J2) as expected from a random array of 
Potts model. 

The specific heat cy per site can be obtained from the energy u as follows, 

where (u^) and (u)^ are the average and squared average over MC steps for u"^ and u 
respectively. The entropy per site, s(T), is obtained from integration over the specific heat 
where s(Ti) and s{T2) are the entropy at the lower and higher temperatures respectively. 
We chose T2 in the high temperature limit where we can take s(T2) = kB^ogeS and used the 
simulation to evaluate the entropy in the ground state, s(0). 

- s{T,) = £^ 'fdT. (6) 

The errors in the determination of the ground state entropy are estimated as follow: (i) by 
comparing the simulated values with the exact result obtained from the order parameter, 
(ii) for the antiferromagnetic phase we compare with published data, (iii) finding the change 



in AS* that result is from changing the parameters in the simulation such as the temperature 
steps, AT, the MC steps or the lattice linear size, L. 




FIG. 3: (a) J1-J2 phase diagram of the orbital structures in simple cubic lattice according to the 
new exchange interaction types Ji and J2 shown in Fig. [2 On the abscissa there is J2, while the 
ordinate represents Ji. (b) The J1-J2 phase diagram, to simplify its study, is divided to region-1 
from the line of the point a = (Ji, J2) = ( — 1, 1) to the line of the point h = (1, 1), region-2 from h 
to c = (1,-1), region-3 from ctod=(— 1,— 1) and region-4 from d to a. For more details see the 
text. 

III. RESULTS AND DISCUSSION FOR THE REGIONAL PHASES 

The simulations show that there is a different and unique order parameter for the six 
regions as shown in J1-J2 phase diagram in Fig. E^. We describe each of these phases in 
turn. The ground state energy and transition temperature are obtained by both simulations 
and analytic reasoning along the four lines, a—^b,b—>-c,c~*d and d —>■ a as presented 
below. 

A. Phase- 1 

This is the well studied ferromagnetic phase that has an entirely ordered ground state 
configuration. Fig. H^. This phase has a three dimensional order parameter corresponding 
to ordering of x, y or z orbits. The simulations show that the region of stability extends 



7 







X 


X 







(a) 




X y X y 




^ X y X 
X ly X ^ 








X 







(b) 




(e) 




FIG. 4: Schematic illustration for the ground state configurations for the regional phases obtained 
by the MC simulation of our model, (a) Phase-1, entirely ordered in each direction, (b) phase-2 is 
inserted between Li and L2 in the J1-J2 phase diagram. The x- and y-layers are alternating in the 
z direction, (c) phase-3, in the x-y plane, the z-state is ordering antiferromagnetically with y-state 
in X direction and with x-state in y direction, (d) phase-4 the well known AF Potts ground state, (e) 
phase-5 corresponds to the orbital ordering in manganites, (f) phase-6 is alternate ferromagnetic 
sheets and checkerboard layers. In all cases only one configuration is shown. 

from Lq where J2 > and Ji > — J2 to Li where Ji > and J2 > (see figure Efe')- The 
analytic expression for the ground state energy is given by, 



-Ji — 2J2 



(7) 



when Ji > 0, there is no frustration in this phase. When Ji < the optimal energy would 
be Uopt = — 2J2 and so in this case the ground state energy in the ferromagnetic phase is 
not optimal. The stabilization energy of this phase is obtained as, Aui = Mi(oo) — ■u(O) = 
|(Ji + 2J2). Since the ground state is completely ordered, as shown in Fig. E^, we expect 
that the value of the ground state entropy, si(0), vanishes and this was confirmed by our 
simulations with an error equal to ±0.005. In Figs. and [7^ we show the excellent 
agreement between the simulated (solid circles) and analytic (solid line) values for the ground 
state energy, 'Ui(O), as a function of | J1I/J2 along the line a-b and \J2\/Ji along the line b-c, 
where the lines were defined in Fie. Eb- 
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FIG. 5: Ji-dependence of (a) the ground state energy, n(0), and (b) transition temperature, Tc, 
from (Ji, J2) = ( — 1, 1) to (1, 1) along the line a-b. 



Figures. and [7|d show the transition temperature, T^, as a function of J1/J2 along 
the line a-b and J2/ Ji along the line b-c in Fig. IHb- We note that the FM Potts model for 
Ji = J2 is well studied and that value of Tc we obtained from the simulation, Tc = 1.8±0.02 
ks/Ji is in agreement with the value obtained in Ref.— . This point is marked by a circle on 
Figs. and[7|D. We shall return to the subject of the specific heat in phase-1 in the region 
where Ji < in the section dealing with Lg. 

Summarising, we have reported known results for the ferromagnetic phase where it is clear 
that our simulation confirmed the well known result that the phase transition of FM Potts 
model is first order, as seen in the plot of the temperature variation of the internal energy, 
Mi(T), and specific heat, cy(T), in Fig. IHl in agreement with other simulation studies and 
with the prediction of the e expansion^^ but in contrast to the result of the position-space 
renormalizat ion-group calculations^!. 
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T(K) T(K) 



FIG. 6: Monte Carlo simulation of the T-dependence of (a) ground state energy, ui{0), (b) the 
specific heat showing the transition temperature, Tc = 1.8K for Ji = J2 = 1- 

B. Phase-2 

In phase-2, ordered ferromagnetic layers (OFM) exist between Li and L2 for — Ji/2 < 
J2 < 0. In this region J2 has become weakly AF and Ji is still strongly ferromagnetic (see 
figure OK)- The ground state configuration consists of alternating FM layers, for example 
X,Y,X,Y,...., along the z direction as seen in Fig. |3Jd. In this phase, in an X layer, an x 
site has four nearest neighbours which are also x-state, two along the x direction and two 
along the y direction, and zero energy between any two neighbouring planes. In this case 
the ground state degeneracy is 6 as there are three ways to choose the normals to the plane 
and then a further factor of two to form the antiferromagnetic arrangement. The ground 
state energy, ^2(0), is, 

W2(0) = - Ji - J2. (8) 

The energy required to get this phase ordered is, Au2 = — |(^i+2 J2) + (-^i + <^2) = |(2Ji + ^2)- 
This phase has long range order, hence, its ground state entropy, 52(0), is = 0.0 ±0.006 as is 
confirmed by simulation. Since J2 < this energy is not optimal. The lowest energy would 
be — Ji but this is frustrated. Because the specific heat as function of T obtained from the 
simulation for phase-2 is as sharp as that for phase- 1 whose transition is well known, we 
believe that phase-2 has a first order transition. Figs. and[3> show ^2(0) and Tc for this 
phase. The line phase Li occurs when Ui{0) = ^2(0) where Ji > and J2 = at this line. 
The variation of the transition temperature with J2 is also continuous at J2 = 0. 
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FIG. 7: J2-dependence of (a) the ground state energy, n(0), and (b) transition temperature, Tc, 
from (Ji, J2) = (1, 1) to (1, —1) along the line b-c. 

C. Phase-3 



This phase is located when we move from L2 to L3 where < Ji < — J2/2 (see figure EK). 
This occurs as J2 becomes more strongly antiferromagnetic, and it has more complicated 
ordering. For clarity we showed in Fig. HJ:; only the 3/4 of sites that have long range order. 
These form a 3d network of antiferromagnetic chains whose energy is zero. The sites left 
blank have two x states nearest neighbours in the x direction, two y states nearest neighbours 
in the y direction and two z states in the z direction. The energy of a state on the blank 
site is exactly — 2Ji whichever state occupies this site. We call this phase the 'Cage' phase 
because the energy comes from the sites left blank, as seen in Fig. HJ:;, that are in a cage. 
However, only 1/4 of the sites are blanks. So the ground state energy, ^3(0), for phase-3 is, 



^3(0) = i(-2Ji) = -\h. 



(9) 



Because the occupation of 1/4 of the sites may be chosen randomly, the analytic expression 
for the entropy is, 

sMI'^B = \log,3 ~ 0.27465. (10) 

The energy needed to order this phase is, Au^, = — |(^i + 2J2) + |Ji = |(|^i — 2J2). 

Again, we have a frustrated phase because of the competition between the strong AF 
J2 and FM Ji. If it were possible to arrange the Potts model so that each one had an 
interaction energy — Ji with two nearest neighbours and zero interaction with the nearest 
neighbours in plane, the energy per site would be —Ji. Figs. andEK show that there is 
excellent agreement between the simulated and analytic results of the ground energy, m(0), 
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as a function J2/ Ji along the line b-c and Ji/| J2I along the line c-d. In this case, the ground 
state entropy obtained by the simulation, at the point J2/ Ji = —0.6, is 0.27 ± 0.005 (see 
Fig. IHId), and this value is in fair agreement with the value which is predicted analytically 

Eq. (uni. 




T(K) T(K) 

FIG. 8: MC simulations of the T-dependence of (a) the internal energy, u^{T)^ at the beginning of 
the phase-3 region at J2/J1 = —0.6 and the end of the region at J2/J1 = —5 (b) entropy, sz{T), 
at J2/J1 = -0.6. 

Additionally, when the ground state energy of phase-2 is equal to that for phase-3, ^2(0) = 
^3(0), the line phase L2 is obtained. From Eqs. (jH)) and 0, we get J2 = —\Ji which agrees 
with the phase boundary L2 obtained from the simulation. 

This phase has a first order transition as Ji — > 0, as seen from the plot of Uz{T) with 
T{K) in Fig. |H^. The simulated specific heat is similar for this phase to that for phase-1. 
The simulation is able to give an accurate value of the ground state entropy to the analytic 
value in the beginning of the phase-3 region but not at the end of this phase where the 
transition temperature is approaching zero. 

It is shown that Tc, obtained from the simulation, as function of J2 for phase-3 region in 
Fig. [7|d, increases dramatically and Tc as function of Ji, Fig. Eb, decreases promptly to zero 
at L3 where there is no a phase transition. The value of Tc is obtained from the simulation 
and is shown in Fig. [Tjp for —1 < J2/J1 < —0.5 and in Fig. for —1 < J2/J1 < 0. The 
transition temperature falls to zero as Ji approaches zero at L3. 
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D. Phase-4 



This is the 'AF' phase which is well known and has been studied extensively2ii22iSS. This 
phase occupies the whole quadrant where both Ji and J2 are AF and it is located between 
L3 at the point (0,-1) and L4 at the point (-1,0), see Ji — J2 phase diagram in Fig. 
Each site in this phase is surrounded with different orbitals (see figure Eli), the ground state 
energy through the whole phase is 

M4(0) = 0, (11) 

where the stabilization energy of this phase is, AM4 = —\{Ji + 2J2) + = ~\{Ji + 2J2). 
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FIG. 9: Ji -dependence of (a) the ground state energy, n(0), and (b) transition temperature, Tc, 
at 1 > Ji > — 1 and J2 = —1. 



The order is understood if the lattice is divided into two sublattices, where one of the 
three states is on the first sublattice and the other two states are distributed randomly on 
the second sublatticei^. This leads to a ground state entropy per site of \kBloge2. But, 
sometimes, at T = 0, states of the lattice are on the wrong sublattice if the surrounding 
neighbours permit it^"^. Accurate Monte Carlo estimates including finite size scaling have 
evaluated^^, 54(0) = 0.3673A;b. This is higher than the estimated, 54(0) = ^ksloge^, by an 
amount = 0.0207A;b. We can quantify the argument of Banavar et al (1982)^^ as follow. The 
probability that a site on the ordered sublattice has six identical neighbours is 2(i)^. In this 
case the site on the ordered sublattice can take one of two values. This gives an analysis for 
the entropy, 

54(0) = ^loge2 + ^%e2 ~ 0.3682fcB. (12) 
This is in good agreement with Wang et a^. There are six ways of defining the ordered part 
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of the ground state coming from two ways of defining the sublattice and three choices of the 
orbit that orders. We find the same ground state configurations throughout this phase and 
the value of the ground state entroDV also takes the same value throuehout this phase. 



u I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — r 




T(K) 

FIG. 10: The internal energy, U4{T), vs T{K) at (Ji, J2) = (-0.1, -1), (-1, -1) and (-1,-0.1). 

There is good agreement between the simulated results of the ground state energy and 
the analytic results as function of J1/IJ2I along the line c-d and ^2/1-^^1! along the line d-e, as 
seen in Figs. andlTTk. Fig. El seems to imply that Phase-4 has a continuous transition 
but the possibility of very weak first order transitions cannot be excluded^"^ because the 
simulation was done on a small cluster. 

E. Phases -5 and -6 

The region where Ji is negative and < J2 < 1 is very interesting because it is divided 
into two related phases. The phases have exactly the same ground state energy but very 
different configurations. The ground state configurations are shown in Figs. 0^ and^f. 

In phase-5 we have an x, y checkerboard pattern in the x — y plane and the planes are 
stacked so that the x and y states are in ordered chains up the z axis. This is the pattern 
of orbits seen in LaMnOs ( see Fig. so we call it LMO phase. The contribution to 
the energy comes from these ferromagnetic chains so the ground state energy is given by 

M5(0) = -J2. 

In phase-6 we have an alternation up the z direction of an x-y checkerboard layer with 
a layer that is occupied by z orbits. We call this the FM-CB phase. The contribution 
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to the ground state energy comes from the ferromagnetic layers. Each z orbit has four z 
nearest neighbours in the x and y directions and so the energy per site of the phase is — 2 J2- 
However, only half of the planes are ferromagnetic so the total energy per site is half the 
energy from the planes which is the same as phase-5. The analytic ground state energy per 
site in phase-5 and phase-6 and the stabilization energies are given below, 

uM = «6(0) = -J2. (13) 

Am5(0) = Au,{0) = -i(Ji - J2). (14) 

Both of these are frustrated in the sense that the lowest possible ground state energy of — 2 J2 
is not accessible. The Monte Carlo calculations give the ground state energy in agreement 
with Eq. (fTSj) as shown in Fig. ITTk . 

The region occupied by phase-5 needs either more Monte Carlo steps or a larger lattice 
size in order to reach a pure LMO ground state. If the system does not come into equilibrium 
then a mixed ground state occurs that includes some of phase-6. 




FIG. 11: (a) J2-dependence of the ground state energy, m(0), at — 1 < J2 < 1 and Ji = —1, (b) 
T-dependence of the internal energy, u^{T), at the point J2/J1 = —0.4 in phase-5 region. 

Because the ground state energies of phase-5 and phase-6 are equal we need a more 
sophisticated argument to obtain the phase diagram. We compare the free energy of the 
two phases in the limit as the temperature approaches zero. This is an example of a phase 
stabilized by disorder, 'order by disorder ',i^>^ which has been used to study frustrated 
Heisenberg models. We note that we are considering a broken discrete symmetry compared 
with the continuous symmetry problems discussed earlier. At T 7^ 0, we calculate the free 
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energy for each phase individually. We write F{T) = —J2 + AF{T) and evaluate AF. 

AF = -keTlogeZ. (15) 

The partition function, Z, is evaluated from the energies of excitations away from the 
ground state where [3 = l/fc^T. 

Z = l + ^e-^^\ (16) 

i 

In the case of phase-5 the LMO phase, we can flip an x site to either y or 2; as shown in 
Figs. IT^ and 112b. The changes in energy associated with the ringed site in cases (a) and 
(b) are given, 

Aea = -Ji -J2 + J2 = -Ji, (17) 
Aefe = +J2. (18) 
In both cases there is a factor of J2 coming from the change in energy per site associated 



X y X ^ ^^^^^^^L * 

F y ® y X ^M^W TG)y 

IjJ^x X y X 'y^ ^jJ^'' 

. I 



y X 
X y X 



7^ 



(a) (b) 



FIG. 12: A single defect in the LMO phase. The site marked with a circle has flipped from x (a) 
to y and (b) to z. 

with the bonds in the z direction. We note that Ji is negative in this region, so, both 
energies Ae^ and Ae^ are positive. The energies of flipping a y site to either x or z are the 
same leading to the following expression for Z^, 

Z5 = 1 + 26-^'^ + 2e-^^^. (19) 

This leads to an expression for AF5 in the low temperature limit, 

AF5 = -NkBTloge[l + 2e-^'^ + 2e--^^^] 

^ -2NkBT[e-^''^ + (20) 
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We now consider the FM-CB layer phase. In this case we have four possibihties. We can 
flip a site in the ferromagnetic plane to either x or y as shown in Figs. IT^ and ()13b ) or we 
can flip a site in the checkerboard plane as shown in Figs. IT^ and IT^ . The change in the 
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FIG. 13: A single defect in the FM-CB phase. The site marked with a circle in FM layer has 
flipped from x (a) to y and (b) to z. And, the site marked with circle in CB layer has flipped from 
z (c) to X and (d) to y. 

energy per site associated with the ringed site in cases (a) and (b) is 



Ae„ = - (-2J2) = 2 J2 
Aeb = -J2 - (-2J2) = ^2, 



(21) 



where, in the FM layers, the ground state energy per site is Ua^i, = — 2J2. In cases (c) and 
(d) Ae is 



Aec = - Ji - = -Ji 
A6rf = -(Ji + J2)-0 



[J1 + J2), 



(22) 



where the ground state in CB layer case is zero. We note that in phase-6, Ji < and 
•J2 < \ Ji\, so, all the energies Aea, .., Aed are positive. The partition function for phase-6, 
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^6, is 



Z, = l + e~^'^^ + e-'^^ + e'^^ + e^-'-^'^^^ . (23) 

The free energy AFg is 

^ -NkBT[e~'^-^^^ + e~-^''^ + e^'^ + e^-^'^-^^^^]. (24) 

The stable phase will be the one with the lower free energy. 

AF5 - AFg = -NkBT[2e''''^ + 26"^^^ - e'^-^^^ - 6"^^^ - e"^^^ - e'^^'+-^^^'^] 

= NkBT[{e^^^ - l){e-^'^ - (25) 

Phase-5 and phase-6 exist in a region where J2 > and Ji < 0. The condition for phase-6 
to be stable is that (AF5 — AFq) > 0. This is given by (2 J2 + Ji) > 0. Thus, we find that 
phase-5 is stable for Ji < and < J2 < — Ji/2. The boundary between phases-5 and 
-6 comes at Ji = — 2J2 and phase-6 is stable for — 2J2 < Ji < — ^2- The ground state 
configuration obtained by the MC simulation confirms this analytic result, see Fig. for 
phase-5 and Fig. E]F for phase-6. The ground state of these phases are expected to have zero 
entropy and this is also confirmed by the simulation with an error equal to ±0.02. 

It is clear from the Fig. ITTk that the analytic ground state energy, u{0), as a function 
of J2/\Ji\ agrees very well with the simulated results. In addition, it decreases continuously 
and steadily from L4 to Lq with increasing J2. One can easily notice that the condition 
for L4 which is U4 = M5 is for J2 = and when M5 = uq the line phase L5 occurs at 
(Ji, J2) = (—1, 0.5). When uq = ui the line phase Lq, J2 = —Ji, is found to separate phase- 
6 and phase-1. However, the transition temperature, Tc, for phase-5 increases slightly, as 
seen in Fig. IT^ . from zero at L4 to 0.29 ± 0.01 at L5, with increasing J2. Fig. ITTb shows 
the T-dependence of the internal energy, u^iT), and Fig. [TUd shows the entropy, s^iT), at 
the point (Ji, J2) = (— 1,0.4). We compare the behaviour of specific heat for phase-5 and 
the upper transition for phase-6 with that for phase-1 to deduce that both have a first order 
transition. 

The phase diagram in phase-6 region is highly unusual for a Potts model because there 
are two district ordering temperatures. These are seen clearly in the specific heat as shown in 
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FIG. 14: (a) J2-dependence of the transition temperature, Tc, from (Ji, J2) = (—1, —1) to (—1, 1), 
(b) T-dependence of the entropy, S5(r), at the point (Ji, J2) = (—1,0.4) in phase-5 region. 

Fig. ITHb . We looked at configurations from the Monte Carlo simulations in the intermediate 
phase in phase-6 region. These had long range order developing in alternate ferromagnetic 
planes, say, in orbital z as shown in Fig. and the intermediate layers were disordered but 
contained predominantly orbitals x and y. 

This may be understood analytically. The ground state energy for the FM layers is 
mfmIO) = — 2J2, the energy needed to order this layer is 

I^ufm = ufm{0) - u{oo) 

= -2J2 + ^(-^i + 2^2) = ^(^1 - 4^2), (26) 

The ground state energy for the CB layers only is zero, so, the energy required to these 
layers to be ordered is 

AucB = mcb(0)-m(oo) 

= + ^(Ji + 2J2) = ^(-^1 + 2^2). (27) 

We check the difference between, Aucb and Aufm as follows, 

AucB - Aufm = i(Ji + 2J2) - ^{Ji - 4J2) = 2J2. (28) 

Since J2 > 0, Aucb > Aufm- This means that the FM layers order first at Tci, before the 
CB layers order at Tc2, where Tci > Tc2 as seen in Fig. IT^ . 

The unusual feature of this phase is that when the ferromagnetic planes are formed the 
orbits on the intermediate planes are free to order independently. The x and y orbitals on 
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T(K) T(K) 

FIG. 15: T-dependence of (a) the internal energy, ue{T), and (b) the specific heat, cv(T), at the 
point (Ji, J2) = (—1,0.7) in phase-6 region. 

the sites on the intervening planes order antiferromagnetically at Tc2- We expected that this 
transition should belong to the class of two dimensional Ising models and should be second 
order. We believe that our simulations are in agreement with this conjecture because the 
peak in the simulated T-dependence of the specific heat for phase-6 at the lower transition 
temperature (T^g), see Fig|15b. in less pronounced than that for the first order transitions we 
have studied. The peak is similar to that for the antiferromagnetic phase (phase-4) which is 
obtained to be second order . There is no reason for the antiferromagnetic order parameter 
to be coherent up the z direction so this phase would show disorder scattering down to low 
temperatures. However, the entropy per site would vary as ksL^'^loge^ and would vanish in 
the thermodynamic limit. This is confirmed by the simulations. 

IV. RESULTS AND DISCUSSION FOR THE BOUNDARY LINES 

The thermodynamics of the boundary lines between two phases will differ from the 
phases on each side because in this case two different configurations can occur with the 
same ground state energy. Because of the extra allowed configurations boundary lines, 
the entropy will be greater or equal to that of both the adjoining phases in all cases. We 
investigate each boundary phase. 
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FIG. 16: Schematic shape for the hne phases ground state ordering (a) 'RFM' layers phase along 
Li which is inserted between phase-1 and phase-2, (b) 'Wood Pile' phase along L2 which is inserted 
between Phase-2 and phase-3, (c) FM and CB layers along L5 alternate randomly with each other 
in X direction, so it is called 'RFMCB' phase, (d) 'FM-disorder' layers phase along Lq which is 
located between phase-6 and phase-1. 

A. Phase Li 

The phase on the line Li in the Ji — J2 phase diagram is investigated. This line separates 
phase-1 ('FM' phase) and phase-2 ('OEM' Layers) where Ji > and J2 = 0. This phase is 
a transition from ferromagnetically coupled layers to antiferromagnetically coupled layers. 
In this case we expect the planar order to be preserved on this line. Fig. ITBk shows that Li 
consists of different (x-layer and y-layer) FM layers alternating randomly with each other 
in z direction. We call this phase 'Random FM' Layers or 'RFM'. Each site has, in-plane, 
two nearest neighbours with Ji interaction and two nearest neighbours with J2 interaction. 
Then, the analytic formula representing the ground state energy per site for this phase for 
J2 = is 

ulM = -Ji- (29) 
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Its ordering energy is, Am^^ = — Ji + |(Ji + 2J2) = |(J2 — Ji). Because the Li phase is 
completely ordered in two dimensions, its ground state entropy per site, Sij(O), tends to 
vary as L~^loge2 and is zero in the thermodynamic limit. The simulated value is ± 0.008. 

It is clear from Fig. that the simulated energy agrees with the analytic value for this 
line phase. The value of the transition temperature is continuous across Li. As seen from 
the simulated plot of T-dependence of specific heat per site, cyiT), in Fig. [T7k the Monte 
Carlo simulations indicate that the transition at Li is similar to that in the ferromagnetic 
phase and hence is expected to be first order. 



B. Phase L2 

This line at Jij J\ = —1/2 separates the antiferro magnetic layer phase from the cage 
phase. On this line the AF effects of J2 prevent the formation of the FM layers and the 
competition with the FM effect of Ji yields FM chains perpendicular to AF chains as seen 
in Fig. ITHb . 





"0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 
T(K) 



FIG. 17: T-dependence of specific heat, cy(T), (a) for the line phase Li, (b) for the line phase L2. 

This phase is obtained when ^2(0) = ^3(0). We get this phase when J2 = —jJi, any 
site could have two nearest neighbours with energy — J2 per site and other two nearest 
neighbours with energy — Ji per site. The total ground state energy is, ul2(0) = ~Ji — J2, 
but we know that J2 = —^Ji, hence. 



^^2(0) = -^Ju 



(30) 



and its ordering energy is, Am^j = ~ki-Ji + 2-^2) + k-Ji = —^(2^2 — k-Ji)- 



22 



This phase is called wood pile as there are ferromagnetic chains of x, y and z orbitals 
running along the x, y and z directions respectively. This leaves one quarter of the sites 
disordered so the ground state entropy is given by 5^3(0) = ^loggS ~ 0.2746A;b which is 
confirmed by the simulations with an error equal to ±0.001. 

The simulated and analytic ground state energy ul2(0) shown in Fig. agree. The 
transition temperature is continuous across L2, Fig. [7|d. The plot of the specific heat shown 
in Fig. [T7b indicates a first order transition compared with the transition of phase-1. 

C. Phase L3 

This phase is at Ji = and J2 < and separates phase-3 and phase-4. The energy of 
this phase is found from this condition ^3(0) = ^4(0) = 0. This is confirmed by simulation 
as shown in Fig. IHK- There is no long range order hence no transition on this line which 
separates two phases having dissimilar long range order, phase-3 and phase-4. However, 
very short FM chains can be seen diagonally with short AF chains along the three axes. 
These different ways of ordering make the ground state entropy higher than that for long 
range ordered phases adjacent to this line. 

While, there is no analytic expression for the ground state energy and entropy, the 
simulated value for the ground state energy, ul.^{0) is equal to that for phase-4, ^4(0), 
in Eq. (Illj) as seen in Fig. The entropy value obtained by the simulation is, 

5^3(0) = (0.423 ± 0.003)/cb. Fig. shows that there is discontinuity in the Ji-dependence 
of Tc at L3. In the phase-3 region, Tc decreases with decreasing of Ji/| J2I and goes to zero 
at L3, but suddenly jumps to a finite value in phase-4, 'AF' phase, and increases linearly 
with decreasing of Ji/\J2\ to become equal to 1.2 A;_b/| Ji| at Ji = J2. 

D. Phase L4 

The phase line at Ji < and J2 = is located between phase-4 and phase-5 (see Fig. 
Et)- The energy on this line, ul^{0) = 0, is confirmed by the simulations as shown in Fig. 
ITTk . There is no broken symmetry in this phase and no evidence of a transition from the 
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specific heat simulations. 

This phase has short range order such that no head-to-head orbital pairs occur. This gives 
a good account of the orbital ordering in a LaMnO^ crystal above its phase transition^. 
The entropy in this phase takes the highest value in the whole phase diagram, ^^^(O) = 
(0.590 ± 0.002)A;5, which is comparable with that obtained experimentally for LaMnO^ by 
Sanchez et af ^. In this case the plot of Tc as a function of J2I J\ is continuous as shown in 
Fig. IT^ and takes its minimum value Tc = on L4. 



FIG. 18: T-dependence of specific heat, cy(T), for the line phase L5. It seems to have a strongly 
first order transition. 

E. Phase L5 

This is a triple point where three transition lines meet as shown in Fig. IT^ and is found 
to be strongly first order as seen from the divergence of the specific heat in Fig. ITHl compared 
with that for phase-1. The ground state energy of L5 is the same as for phase-5 and phase-6. 



This phase is a mixed between the two distinct orderings of phase-5 and phase-6. As these 
phases both have entropies that tend to zero in the thermodynamic limit the entropy is 
expected to vanish on L5 too. This is confirmed by the simulation with an error equal 
to ±0.006. The transition temperature is shown on Fig. ^1 It is seen that it lies on a 
continuous of the lines from phase-5 and phase-6. 




T(K) 



Mi,(0)=M5(0)=M6(0) = -J2. 



(31) 
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F. Phase Le 



This line separates phase-6 (FM-CB) and phase-1 (FM) phase and occurs for uq{0) = 
Ui{0) which is found from Eq. ((7j) and (jl3j) to occur for Ji + J2 = 0. Phase-6 has two transi- 
tion temperatures. The lower line, Tc2, goes to zero on Lq. The ground state configuration 
of the phase on Lq, as shown in Fig. ITBH . is ferromagnetic layers alternating with disordered 
layers. The phase of Lg is called 'FM-Disorder' phase. 




T(K) 

FIG. 19: T-dependence of specific heat, cy(T), at Ji = —1, —0.9, —0.8, and -0.7 and J2 = 1.0. The 
upper three panels correspond to phase-1 and the bottom panel is for Lq. There exists a short 
range order peak joining with the transition peak between Lq and Ji = —0.7. 

The ground state energy can be obtained analytically. The disordered layers have total 
energy equal to zero but the FM layers are ordered with, say y states in x-z plane. So the 
exchange interaction for each site in a FM layer is — 2J2. The total energy at the ground 
state for the whole lattice is ulq{0) = — J2- We divided by 2 because only the half of the 
layers are FM. However, the ground state entropy of the phase on this line comes from the 
disordered layers which have two states distributed randomly. We predict Sig(O) = ^loge2, 
again, we divided by 2 because the entropy comes only from the disordered layers which 
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comprise half of the lattice. The simulated value of the ground state entropy agrees with 
the analytic one with an error equal to ±0.005. 

We consider phase-6, phase-1 and the phase on the line joining them, Lg, together. The 
specific heat at the upper transition evolves from a sharp anomaly as shown in Fig. IT^ to 
a broader peak at Lg as shown in Fig. E| In phase-1 the peak appears to correspond to the 
onset short range order. 

The fact that we appear to have a line of transition at Tc^ in phase-6 that evolves into 
a short range order peak in the FM phase (phase-1) is reminiscent of the critical end point 
seen in the liquid-gas transition. However, we see none of the critical phenomena associated 
with a second order transition. A finite system is unable to distinguish between short range 
order and long range order below the upper critical temperature on Lg. This area of the 
phase diagram needs further investigations. 



TABLE I: The region and line phases, the analytic ground state energy Un{0) and entropy Sn(0) 
per site and transition temperatures Tc obtained with Sn(0) by MC simulation of three states 3d 
AAFP model for cubic lattice with L = 8. 



Phase No. 


Phase Name 


analy. tin(O) 


analy. s„(0)/fcB 


simulat. Sn{0)/kB 




phase-1 


FM 


- Jl - 2J2 





0.0 ±0.005 




phase-2 


OFM 


- Jl - J2 





0.0 ±0.006 




phase-3 


Cage 


-Ji/2 




0.270 ± 0.005 




phase-4 


AF 


0.0 


\loge2 + ^loge2 


0.3673'' 




phase-5 


LMO 


-J2 





0.0 ±0.02 




phase-6 


FM-CB 


-J2 





0.0 ±0.02 




Li (J2 = 0) 


RFM 


-Jl 





0.0 ±0.008 


0.48Ji 


L2 (J2 = -Ji/2) 


Wood Pile 


-Ji/2 


lloge3 


0.274 ± 0.001 


O.39J1 


L3 {Jl = 0) 


Disorder 







0.423 ± 0.003 




L4 ( J2 = 0) 


AAFP 







0.590 ± 0.002 




L5 (Jl = -J2/2) 


RFMCB 


-J2 





0.0 ± 0.006 


O.29J2 


Lg (Jl = -J2) 


FM-Disorder 


-J2 


lloge2 


0.342 ± 0.005 


O.57J2 



^Ref {^) 
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All six regional and line phases obtained pertaining to this phase diagram have been 
analysed. Two of the phases were the well studied ferromagnetic and antiferromagnetic 
Potts models with Ji — J2 — J and J > and J < respectively and our results confirmed 
the known expressions in these cases. One of the phases corresponds to the ordering seen in 
LaMnO^. 

We used a combination of Monte Carlo simulations and analytic reasoning to obtain our 
results. The Monte Carlo simulations were run from high to low temperatures and the 
configurations obtained at the lowest temperatures were analyzed to find if the symmetry 
had been broken and if so to identify the order parameter. This identified the existence 
of a phase transition. The ground state energy was obtained from the simulations for all 
phases. We used the observed configuration to recalculate the ground state energy and in 
all cases obtained excellent agreement with the simulated values. The analytic expressions 
for the ground state energies also enabled us to identify the stability lines for each phase by 
equating the energies of two neighbouring phases. In all cases the simulations confirmed the 
phase stability lines that we had found analytically. In most cases knowledge of the ground 
state configuration enabled us to obtain an analytic expression for the entropy in the ground 
state. This was harder to obtain from the simulations particularly, the simulations indicated 
that there was a first order phase transition but again the results of the simulations agreed 
with the analytic results to within the errors. In the case of two boundary lines there was 
no ordering and so there was no analytic expression for the entropy and the only estimate 
was obtained from the simulations. 

The transition temperatures were obtained from the simulations and plots were presented 
of the variation of the transition temperatures with the variation of the coupling constants. 

It was found that only the ferromagnetic and antiferromagnetic Potts phases were not 
frustrated in all other cases the ground state energy was higher than the optimal value. A 
number of novel phases are obtained from our phase diagram. The most unusual results are 
listed below. 

Phase-3 (or 'Cage' phase) has the most unusual ground state configuration in the phase 
diagram. Three quarters of the sites are ordered but a large contribution to the ground 
state stabilization energy comes from the disordered sites. It is seen that the transition 
temperature falls to zero as Ji is reduced. The energy stabilization from the disordered sites 
varies linearly with Ji. 
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Phase-5 and phase-6 have the same ground state energy but with very different config- 
urations. Their regions of stabihty were found from considering the free energy at finite 
temperature which was evaluated using a calculation of the elementary excitations. This 
was an example of 'order by disorder' that had previously been applied to models with 
continuous symmetr y^^'^^ . The method gave the ordering of the phases correctly and also 
identified the line between them. 

Phase-6 has two well-defined phase transitions. This is very unusual for a Potts model. 
These phases were found in the simulations and an analytic discussion was given that added 
understanding. The temperature of the two transitions coincided at the boundary with 
phase-5 giving a triple point. At the other end of the phase stability the lower transition 
temperature went to zero at the boundary with the ferromagnetic phase and the upper 
transition evolved into the onset of short range order. 

In some cases the phase line separated two ordered phases that had various elements 
in common, an example of that was the transition between ferromagnetism layers stacked 
antiferromagnetically. In this case the common feature, ordered layers, was preserved on 
the boundary line and the transition temperature on the boundary line stayed finite and 
the ground state entropy was zero. In other cases the adjoining phases had no elements in 
common and in these cases the there was no transition at any finite temperature on the 
boundary line and a high value for the entropy was found in the low temperature limit. 

In summary this simple Potts model shows a great diversity of phases and different critical 
behaviour. At least one phase corresponds to a physically realized orbital ordering and it 
will be interesting to see if any of the other more exotic phases have a physical realization. 
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Corrections 

• The abstract has been changed shghtly. 

• Page 2: The subtitle 'Potts model' in the introduction has been omit- 
ted. 

• Page 4 and 5: We apologise for the typo in the Hamiltonians given in 
equations 2 and 3. These have rewritten using the notation of Wang 
(ref. 28). 

• Page 25: Phase-1 close to Lg has been presented at greater length. 
Figure 19 has been added. 

• Table 1 has been moved through the conclusion section with some mod- 
ification. 

• In the L5 subsection, we changed the phrase "tricritical point" to be 
"triple point". 

• In the phase-3 subsection, the sentence "the finite order transition tem- 
perature" has been changes to be "the finite transition temperature". 

• We corrected "For clarify" to be "For clarity" and changed "These are 
from" to be "These form" in phase-3 subsection. Finally, we corrected 
the caption of figure 6. 
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